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Abstract. In this talk I review NRQCD predictions for the production of charmonium at the 
Tevatron. After a quick presentation of the NRQCD factorization formalism for production and 
decay I review some old results and discuss how they compare to recent data. Following this I 
discuss some recent work done with Adam Leibovich and Ira Rothstein. 


INTRODUCTION 

Heavy quarkonia have proven fruitful in helping us gain a better understanding of QCD. 
Early theoretieal analyses of quarkonia decay were based on the color-singlet model 
(CSM) [1]. The underlying assumption of this model is that the heavy-quark-antiquark 
pair has the same quantum numbers as the quarkonium meson. (For example the bb 
that forms an T must be in a color-singlet configuration.) One consequence of such 
a restrictive assumption is that theoretical predictions based on the CSM are simple, 
depending on only one nonperturbative parameter. However, the CSM does not provide 
a systematic approach to studying quarkonium. This is clear in E-wave decays where 
infrared divergences signal the breakdown of the CSM [2]. 

In order to systematically study nonrelativistic systems long distance physics needs 
to be separated from short distance physics. This can be accomplished with a proper 
effective field theory which provides a power counting that determines relevant opera¬ 
tors. In most effective theories the power counting is based upon dimensional analysis, 
however, for non-relativistic QCD (NRQCD) [3, 4] this is not the case. Instead, it is an 
expansion in the parameter v, the relative velocity of the heavy quarks. This leads to 
the result that operators of the same dimension may be of different orders in the power 
counting. Inclusive decay rates and production cross sections are now understood within 
the framework of NRQCD factorization [3], where decay rates and production cross sec¬ 
tions are predicted in a systematic double expansion in and v. These predictions have 
met with varying degrees of success. 

In the first half of my talk I give a quick review of NRQCD, and the NRQCD 
factorization formalism for production and decay. This gives the background needed 
to understand theoretical predictions for production and decay of J/\\r and \|/'. I do not 
attempt to review all of the predictions, instead I focus on the transverse momentum 
distribution of J/\\f and t|/' at the Tevatron. For the reason that here lies the earliest 
success and the greatest challenge to NRQCD factorization. In Ref. [5] the ‘anomaly’ 
(a factor of 30 discrepancy between the CSM prediction and date) was resolved using 
NRQCD. However, the initial data on the polarization of these states at large transverse 
momentum [6] seems to be at odds with the NRQCD prediction. 
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In the second half of the talk I will discuss a possible resolution to the charmonium 
polarization puzzle at the Tevatron. This is based on recent work [7] with Ira Rothstein 
and Adam Leibovich wherein we propose an alternative power counting for charmo¬ 
nium. I do not present all the pieces of evidence which seems to tell us that the effective 
field theory which best describes the 7/v[t system may not be the same theory which best 
describes the T. For that the reader is directed to the literature. I merely give a quick 
review of the new power counting, and then proceed to discuss how this changes the 
predictions for 7/\|/ and \|/' polarization at the Tevatron. 


NRQCD 


The power counting depends upon the relative size of the four scales (m, mv, mv^, Aqcd) • 
If we take m > mv > mv^ ~ Aqcd the bound state dynamics will be dominated by ex¬ 
change of Coulombic gluons with (E ~ mv^^p = mv). This hierarchy has been assumed 
in the NRQCD calculation of production and decay rates and is probably a reasonable 
choice for the T system, where mv ~ 1.5 GeV. However, whether or not it is correct for 
the J /\\f, where mv 100 MeV remains to be seen. 

The power counting can be established in a myriad of different ways. Here I will 
follow the construction of [4], which I now briefly review. There are three relevant 
gluonic modes [8]: the Coulombic (mv^,mv), soft (mv, mv) and ultrasoft (mv^,mv^). 
The soft and Coulombic modes can be integrated out leaving only ultrasoft propagating 
gluons. In the process of integrating out these modes large momenta must be removed 
from the quark field. This is accomplished by rescaling the heavy quark fields by a factor 
of exp (ip ■ x) and labeling them by their three momentum p. The ultrasoft gluon can only 
change residual momenta and not labels on fields. This is analogous to HQET, where the 
four-velocity labels the fields and the nonperturbative gluons only change the residual 
momenta [9]. This rescaling must also be done for soft gluon fields [10] which, while 
they cannot show up in external states, do show up in the Lagrangian. After this rescaling 
a matching calculation leads to the following tree level Lagrangian [4] 


L = 
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where we have retained the lowest order terms in each sector of the theory. The matrices 
and are the color matrices for the 3 and 3 representations, respectively. Note the 
last term is the Coulomb potential, which is leading order and must be resummed in the 
four-quark sector, while the other non-local interactions arise from soft gluon scattering. 

All the operators in the Lagrangian have a definite scaling in v, and the spin symmetry, 
which will play such a crucial role in the polarization predictions, is manifest. The two 
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subleading interactions which will dominate my discussion are the electric dipole (El) 

=¥i--AVp, (2) 

and the magnetic dipole (Ml) 


ft 
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^M1 = • 

The El interaction is down by a factor of v while the Ml is down by a factor of v^. 


( 3 ) 


NRQCD FACTORIZATION FORMALISM 

In the NRQCD factorization formalism developed in Ref. [3] decay rates and production 
cross sections are written as a sum of products of Wilson coefficients encoding short 
distance physics and NRQCD matrix elements describing long distance physics. In this 
formalism a general decay process is written as 

Ty/v = I I V). (4) 

The matrix element represents the long distance part of the rate and may be thought of 
as the probability of finding the heavy quarks in the relative state n, while the coefficient 
Cis+\ij{m^<Xs) is a short distance quantity calculable in perturbation theory. The sum 
over operators may be truncated as an expansion in the relative velocity v. Similarly, 
production cross sections may be written as 

d<5 = I O" I 0). (5) 

n 

Here the short distance cross section for a reaction involving two 

partons, i and j, m the initial state, and two heavy quarks in a final state, labeled by 
n, plus X. This part of the process is calculable in perturbation theory, up to possible 
structure functions in the initial state. The production matrix elements, which differ from 
those used in the decay processes, describe the probability of the short distance pair in 
the state n to hadronize, inclusively, into the state of interest. The relative size of the 
matrix elements in the sum are again fixed by the power counting which we will discuss 
in more detail below. 

The formalism for decays is on the same footing as the operator product expansion 
(OPE) for non-leptonic decays of heavy quarks, while the production formalism assumes 
factorization, which is only proven, and in some applications of production this is not 
even the case, in perturbation theory ^ [11]. The trustworthiness of factorization depends 
upon the particular application. I have reviewed these results here to emphasize the point 
that when the theory is tested one is really testing both the factorization hypothesis as 
well the validity of the effective theory as applied to the J/\\f system. Thus, care must be 
taken in assigning blame when theoretical predictions do not agree with data. 


^ For a discussion of factorization in NRQCD see Refs. [3, 12]. 
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TABLE 1. Scaling of matrix elements 
NRQCDfo and NRQCDc. 

relevant for \|/ production in 


(0T(35i)) 


(^?I('5o)) 


NRQCDfo 





NRQCDc 

(AgcD/wie)** 

(Aqcd/»Jc)"^ 

(Aqcd/»Jc)^ 

(Aqcd/wJc)"^ 


CHARMONIUM PRODUCTION AT THE TEVATRON 

Having introduced NRQCD and the factorization formalism I will now turn to theoretical 
predictions of 7/\|/ and \\r' production at the Tevatron. The leading order in v contribution 
to J/\\f production is through the color-singlet matrix element (6)^(^5i)), since the 
quantum numbers of the short distance quark pair matches those of the final state. All 
other matrix elements need insertions of operators into time ordered products to give 
a non-zero result, and are therefore suppressed compared to the color-singlet matrix 
element above. For instance, the matrix element {Oj{^So)) vanishes at leading order. 
The first non-vanishing contribution comes from the insertion of two Ml operators into 
time ordered products, thus giving a suppression. The scaling of the relevant matrix 
elements for \\f production are shown in Table 1 under NRQCD^ (for reasons which I will 
explain later). It appears from just the v counting that only the color-singlet contribution 
is important. However, other contributions can be enhanced by kinematic factors. At 
large transverse momentum, fragmentation type production dominates [13], and only 
the (6)^(^Si)) contribution is important. Without the color-octet contributions {i.e., the 
Color-Singlet Model), the theory is below experiment by about a factor of 30. By adding 
the color-octet contribution the fit to the data is very good [5]. 

Once the color-octet matrix elements are fit to the unpolarized date it is possible to 
make a parameter free prediction for the polarization of 7/\|/ and V the Tevatron: 
they are predicted to be transversely polarized at large pr- This is because at large 
transverse momentum, the dominant production mechanism is through fragmentation 
from a nearly on shell gluon to the octet state. The quark pair inherits the polarization 
of the fragmenting gluon, and is thus transversely polarized [14]. The leading order 
transition to the final state goes via two El, spin preserving, gluon emissions. Higher 
order perturbative fragmentation contributions [15], fusion diagrams [16, 17], and feed- 
down for the 7/t|/ [18] dilute the polarization some, but the prediction still holds that as 
Pt increases so should the transverse polarization. Indeed, for the \\f', at large pj ^ nic, 
we expect nearly pure transverse polarization. This prediction seems to be at odds with 
the initial data which seems to suggest that the J/\\r and the are unpolarized or slightly 
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longitudinally polarized as pT increases [6].^ If, after the statistics improve, this trend 
continues one is left with two obvious possibilities: 1) The power counting of NRQCD 
does not apply to the 7/\|/ system. 2) Factorization is violated “badly”, meaning that 
there are large power corrections. 


NRQCDc 

In this section I discuss the work done with Adam Leibovich and Ira Rothstein in Ref. [7]. 
In that work we marshaled evidence that the NRQCD power counting might not apply to 
the 7/\|/ system. We did not consider the second possibility mentioned at the end of the 
previous section: that factorization is violated. I will not be considering this possibility 
either. 

The standard NRQCD methodology, which is based upon the hierarchy m > mv > 
mv^ ~ Aqcd, has been applied to the 7/\|/ as well as the T systems. While it seems quite 
reasonable to apply this power counting to the T system, it is not clear that it should 
apply to the J/\\f system. Indeed, I believe that the data is hinting toward the possibility 
that a new power counting is called for in the charmed system. 

If NRQCD does not apply to the J/\\f system, then one must ask: is there another 
effective theory which does correctly describe the 7/\|/? One good reason to believe that 
such a theory does exist is that NRQCD, as formulated, does correctly predict the ratios 
of decay amplitudes for exclusive radiative decays. Using spin symmetry the authors of 
[14] made the following predictions: 

r(XcO ^ j/^ + Y) : r(Xci Y/V + Y) : ^ J/^ + Y) : rje + y) 

= 0.095 : 0.20 : 0.27 : 0.44 (theory) 

= 0.092±0.041: 0.24± 0.04 : 0.27 ±0.03 : unmeasured (experiment). (6) 

Thus, an alternative formulation of NRQCD must preserve these predictions yet yield 
different predictions in other relevant processes. 

Let me now consider the alternate hierarchy m > mv ~ Aqcd- One might be tempted to 
believe that in this case the power counting should be along the lines of HQET, where the 
typical energy and momentum exchanged between the heavy quarks is of order Aqcd- 
However, this leads to an effective theory which does not correctly reproduce the infra¬ 
red physics. With this power counting, the leading order Lagrangian would simply be 

^HQET = VjY)oV[tv, (7) 

where the fields are now labeled by their four velocity. This is a just a theory of time-like 
Wilson lines (static quarks) which does not produce any bound state dynamics. Thus I 
am forced to conclude that the typical momentum is of order Aqcd, whereas the typical 
energy is so that D^/ (2m) is still relevant. I will call this theory NRQCD^, and 

will refer to the traditional power counting as NRQCD^ as I assume that it describes the 
bottom system. 


^ The data still has rather large error bars, so we should withhold judgment until the statistics improves. 
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The power counting of this theory is now along the lines of HQET where the ex¬ 
pansion parameter is Aqcd/^q- However the residual energy of the quarks is order 
'^qcd/^2’ while the residual three momentum is Aqcd- Thus one must be careful in the 
power counting to differentiate between time and spatial derivatives acting on the quark 
fields. As far as the phenomenology is concerned, perhaps the most important distinction 
between the power counting in NRQCDc and NRQCD^, is that the magnetic and electric 
gluon transitions are now of the same order in NRQCDc. This difference in scaling does 
not disturb the successes of the standard NRQCD^ formulation but does seem help in 
some of its shortcomings. 


NRQCDc PREDICTIONS 

The relative size of the different matrix elements change in NRQCD^. In particular, the 
Ml transition is now the same order as the El transition. The new scaling is shown in 
Table 1. ^ Due to the dominance of fragmentation at large transverse momentum, we 
need to include effects up to order /nicY, since the (O^i^Si)) matrix element 

will still dominate at large pj- 

Is this consistent? The size of the matrix elements is a clue. Extraction of the matrix 
elements uses power counting to limit the number of channels to include in the fits. Cal¬ 
culating 7/\|/ and \\f' production up to order (Aqcd/«^c)^ in NRQCD^ requires keeping 
the same matrix elements as in NRQCD^. Previous extractions of the matrix elements 
only involve the linear combination 

M* = {O*('so)) + (8) 

me 

with r ^ 3 — 3.5, since the short-distance rates have similar size and shape. In the new 
power-counting, I can just drop the contribution from {0^(^Pj)), since it is down by 
(Agco/nic)^ ~ I/IO compared to {Oj{^So)). It is the same order as but is 

not kinematically enhanced by fragmentation effects. The extraction from [18] would 
then give for the J/\\r and \|/' matrix elements 

(6>g/^(^5o)) : (6>8/^(^5 i)) = (6.6±0.7) x 10^^ : (3.9±0.7) x 10^^ ^ 17 : 1, 
(OjVSo)) : (6>f(^5i)) = (7.8 ± 3.6) X 10^3 : (3.7 ± 0.9) X 10^3 ^2:1. (9) 

Other extractions have various values of the hierarchy, ranging from 3 : 1 to 20 : 1 [21]. 
While the relation of the color-octet matrix elements in the 7/v[t system is indeed in 
agreement with the NRQCD^ power counting, the does not look to be hierarchical. 
However, it should be noted that the statistical errors in the V|/' extraction, quoted above, 
are quite large. Eurthermore, there are also large uncertainties introduced in the parton 
distribution function. The above ratios used the CTEQ5E parton distribution functions. 
If we take the central values from [18] for the MRST98EO distribution functions, we 


^ These results reproduce those given in [20] when X is taken to be 1 in this reference. 
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find the ratio 3:1. On the other hand, the 7/v[t extraction is much less sensitive to the 
choice of distribution function. Given the statistical and theoretical errors, it clear that 
the t|/' ratio is not terribly illuminating. 

Let me now consider the extraction of these color-octet matrix elements in the T sec¬ 
tor [22], where according to NRQCD^ power counting there is should be no hierarchy: 


00 

('So)) : 


e^l)) 

= (5.4 ± 4.3+^;5) X 10^2 : 

(3.6±1.9+{j) X 10~^ 





1:1, 


(or> 

(‘5o)> : 

00 

e^l)) 

= (-10.8±9.7;|J) X 10 

: (16.4±5.7lJ}) X 10^^ 





^ 1:1, 


(or 

(‘5o)> : 

(Of"' 

e^l)) 

= (13.6 ± 6.81^5^) X 10^ 

2 : (2.0±4.i;|]j) X 10^2 





6:1. 

(10) 


For the T(35) and T(2S) we observe that there is indeed no hierarchy, while for the 
T( 15) it appears there may be a hierarchy. However, it is not possible to draw any strong 
conclusions from these data because the errors on the extractions are large. In fact the 
ratio for the T(15) color-octet matrix elements is 1 : 1 within the one sigma errors. 
Furthermore, these matrix elements are those extracted subtracting out the feed down 
from the higher states. While phenomenologically it is perfectly reasonable to define the 
subtracted matrix elements, I believe that, since the matrix elements are inclusive, one 
should not subtract out the feed down from hadronic decays when checking the power 
counting. In principle this subtraction should not change things by orders of magnitude, 
but nonetheless it can have a significant effect. Indeed, if one compares the ratios for 
inclusive matrix elements, which do not have the accumulated error, then the ratios come 
out to be 1 : 1, even for the T(15) [22]. 

With NRQCDf, the intermediate color-octet ^5i states hadronize through the emission 
of either two E\ or Ml dipole gluons, at the same order in l/nic. Since the magnetic 
gluons do not preserve spin, the polarization of t|/ produced through the {Oj{^Si)) can 
be greatly diluted. The net polarization will depend on the ratio of matrix elements 

Rm/e = ( 11 ) 

lUed^MO I T{Mi{xi)Mi{x2)\\f'^T^aiX)aUHT{Mi{x3)Mi{x4)x^T-aiy\f) \ 0 ) 
lUed^MO I T(El(xl)El(x2)^\f^T^aa)aUHTiEl(x3)El(x4)x'^T^Oi^\f) \ 0 ) 


where 

aJ^a^ = £|// + X)(//-t-X I . (12) 

x 

The leads to the polarization leveling off at large pj at some value which is fixed by 
In Fig. 1, we show the prediction for 7/\|/ and t|/' polarization at the Tevatron. The 
data is from [6]. The three lines correspond to different values for (dashed), 

1 (dotted), oo (solid)). The dashed line is also the prediction for NRQCD/,. The residual 
transverse polarization for 7/\|/ at asymptotically large pj is due to feed down from % 
states. The non-perturbative corrections to our predictions are suppressed by 
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Pt (GeV) 

FIGURE 1. Predicted polarization in NRQCDc for 7/\|/ and \|/' at the Tevatron as a function of pj. The 
three lines correspond to (dashed), 1 (dotted), (solid)). The dashed line is also the prediction 

forNRQCDfo. 


CONCLUSION 

In this talk I have reviewed NRQCD, and the NRQCD faetorization formalism whieh 
is used to make predietions for the produetion and deeay of eharmonium and bottomo- 
nium. I did not diseuss all of these predietions. Instead I foeused on what I believe to 
be the most important predietion of NRQCD faetorization: the transverse momentum 
distribution of unpolarized and polarized J/'^ and \|/' produeed at the Tevatron. Beeause 
the unpolarized data ean be used to determine the unknown eolor-oetet matrix elements, 
it is possible to make a parameter free predietion for polarized produetion. This provides 
a elean test of the NRQCD faetorization formalism. Moreover the quality of the data for 
unpolarized produetion is good, and while the data for polarized produetion has large 
error bars it is expeeted to get better. 

The NRQCD faetorization formalism prediets the 7/\|/ and \|/' to be transversely 
polarized at large pr- This is beeause at large transverse momentum, the dominant 
produetion meehanism is through fragmentation from a nearly on shell gluon to the oetet 
state. The quark pair inherits the polarization of the fragmenting gluon, and is thus 
transversely polarized [14]. The leading order transition to the final state goes via two 
El, spin preserving, gluon emissions. Various eorreetions dilute the polarization some, 
but the predietion still holds that as pr inereases so should the transverse polarization. 
Indeed, for the \\f', at large pj 3> m^, we expeet nearly pure transverse polarization. 
The eurrent experimental results [6] show no or a slight longitudinal polarization, as pj 
inereases. If, after the statisties improve, this trend eontinues, then it will be the smoking 
gun that leads us to eonelude that either NRQCD is not the eorreet effeetive field theory 
for eharmonia, or that faetorization fails in these proeesses. 

The possibility that NRQCD is not the eorreet effeetive theory for eharmonium leads 
me to ask: is there any reason to believe that there is any effeetive theory to eorreetly 
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describe the 7/\|/? I believe that the spin symmetry predictions for the ratio of % decays 
clearly answers this question in the affirmative. Assuming that such an effective theory 
exists, then is it NRQCD^ or NRQCD^? As I have shown the two theories do indeed 
make quite disparate predictions, which in principle should be easy to test. 

However, these tests can be clouded by the issues of factorization and the convergence 
of the perturbative expansion. One would be justified to worry about the breakdown of 
factorization in hadro-production at small transverse momentum. However, for large 
transverse momentum one would expect factorization to hold, with non-factorizable 
corrections suppressed by powers of m^/pr-As far as the perturbative expansion is 
concerned, it seems that for most calculations the next-to-leading order results are 
indeed smaller than the leading order result [23, 24, 25], though, the NNLO calculation 
performed, in the leptonic decay width [26], is not well behaved at this order. 

In the end I believe the data will be the final arbiter. The polarization measurement 
may fall in line with the NRQCD^ prediction. Or the data may result in longitudinal 
polarization for 7/\|/ and \|/', in which case it may be that neither NRQCDe nor NRQCD^ 
are the correct theory. 


ACKNOWLEDGMENTS 

I would like that thank my collaborators Adam Leibovich and Ira Rothstein. Thanks also 
go the organizers of this conference. This work was supported in part by the Department 
of Energy under grant number DOE-ER-40682-143 and DE-AC02-76CH03000. 


REEERENCES 

1. W. . Buchmuller, Amsterdam, Netherlands: North-Holland (1992) 316 p. (Current physics-sources 
and comments, 9). 

2. G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 46, 1914 (1992) [arXiv;hep-lat/9205006]. 

3. G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 51, 1125 (1995) [Erratum-ibid. D 55, 5853 
(1995)1 [hep-ph/9407339]. 

4. M. E. Luke, A. V. Manohar and I. Z. Rothstein, Phys. Rev. D 61, 074025 (2000) [hep-ph/9910209]. 

5. E. Braaten and S. Eleming, Phys. Rev. Lett. 74, 3327 (1995) [arXiv:hep-ph/9411365]. 

6. T. Affolder et al. [CDE Collaboration], Phys. Rev. Lett. 85, 2886 (2000) [arXiv;hep-ex/0004027]. 

7. S. Eleming, I. Z. Rothstein and A. K. Leibovich, Phys. Rev. D 64, 036002 (2001) [arXiv:hep- 
ph/0012062]. 

8. M. Beneke and V. A. Smirnov, Nucl. Phys. B 522, 321 (1998) [arXiv:hep-ph/9711391]. 

9. A. V. Manohar and M. B. Wise, Cambridge Monogr. Part. Phys. Nucl. Phys. Cosmol. 10, 1 (2000). 

10. H. W. Griesshammer, Phys. Rev. D 58, 094027 (1998) [arXiv:hep-ph/9712467]. 

11. See ref. [3] and related proofs in, J. C. Collins, D. E. Soper and G. Sterman, Adv. Ser. Direct. High 
Energy Phys. 5, 1 (1988). 

12. M. Beneke and I. Z. Rothstein, Phys. Rev. D 54, 2005 (1996) [Erratum-ibid. D 54, 7082 (1996)] 
[arXiv:hep-ph/9603400]. 

13. E. Braaten and T. C. Yuan, Phys. Rev. Lett. 71, 1673 (1993) [arXiv:hep-ph/9303205]. 

14. P. L. Cho and M. B. Wise, Phys. Lett. B 346, 129 (1995) [arXiv;hep-ph/9411303]. 

15. M. Beneke and I. Z. Rothstein, Phys. Lett. B 372, 157 (1996) [Erratum-ibid. B 389, 769 (1996)] 
[arXiv;hep-ph/9509375]. 

16. M. Beneke and M. Kramer, Phys. Rev. D 55, 5269 (1997) [arXiv;hep-ph/9611218]. 

17. A. K. Leibovich, Phys. Rev. D 56, 4412 (1997) [arXiv:hep-ph/9610381]. 


9th International Symposium on Heavy Flavor Physics 


February 1, 2008 


9 




18. E. Braaten, B. A. Kniehl and J. Lee, Phys. Rev. D 62, 094005 (2000) [arXiv:hep-ph/9911436]; 
B. A. Kniehl and J. Lee, Phys. Rev. D 62, 114027 (2000) [arXiv:hep-ph/0007292]. 

19. E. Braaten and J. Lee, Nucl. Phys. B 586, 427 (2000) [arXiv;hep-ph/0004228]. 

20. M. Beneke, arXiv;hep-ph/9703429. 

21. P. L. Cho and A. K. Leibovich, Phys. Rev. D 53, 150 (1996) [arXiv:hep-ph/9505329]; ibid. D53, 
6203 (1996) [arXiv:hep-ph/9511315]; 

B. Cano-Coloma and M. A. Sanchis-Lozano, Nucl. Phys. B 508, 753 (1997) [arXiv:hep- 
ph/9706270]; 

C. G. Boyd, A. K. Leibovich and I. Z. Rothstein, Phys. Rev. D 59, 054016 (1999) [arXiv:hep- 
ph/9810364]; 

B. A. Kniehl and G. Kramer, Eur. Phys. J. C 6, 493 (1999) [arXiv;hep-ph/9803256]. 

22. E. Braaten, S. Lleming and A. K. Leibovich, Phys. Rev. D 63, 094006 (2001) [arXiv:hep- 
ph/0008091]. 

23. M. Beneke, I. Z. Rothstein and M. B. Wise, Phys. Lett. B 408, 373 (1997) [arXiv:hep-ph/9705286]. 

24. L. Maltoni and A. Petrelli, Phys. Rev. D 59, 074006 (1999) [arXiv;hep-ph/9806455]. 

25. A. Petrelli, M. Cacciari, M. Greco, L. Maltoni and M. L. Mangano, Nucl. Phys. B 514, 245 (1998) 
[arXiv;hep-ph/9707223]; 

E. Braaten and Y. Q. Chen, Phys. Rev. D 55, 7152 (1997) [arXiv:hep-ph/9701242]. 

26. M. Beneke, A. Signer and V. A. Smirnov, Phys. Rev. Lett. 80, 2535 (1998) [arXiv:hep-ph/9712302]. 

27. E. Braaten and Y. Q. Chen, Phys. Rev. D 57, 4236 (1998) [Erratum-ibid. D 59, 079901 (1998)] 
[arXiv:hep-ph/9710357]; 

G. T. Bodwin and Y. Q. Chen, Phys. Rev. D 60, 054008 (1999) [arXiv:hep-ph/9807492]. 


9th International Symposium on Heavy Flavor Physics 


February 1, 2008 


10 



